Conservative self-organized extremal model for wealth distribution 
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We present a detailed numerical analysis of the modified version of a conservative self-organized 
extremal model introduced by Pianegonda et. al. for the distribution of wealth of the people in a 
society. Here the trading process has been modified by the stochastic bipartite trading rule. More 
specifically in a trade one of the agents is necessarily the one with the globally minimal value of 
wealth, the other one being selected randomly from the neighbors of the first agent. The pair 
of agents then randomly re-shuffle their entire amount of wealth without saving. This model has 
most of the characteristics similar to the self-organized critical Bak-Sneppen model of evolutionary 
dynamics. Numerical estimates of a number of critical exponents indicate this model is likely to 
belong to a new universality class difi'erent from the well known models in the literature. In addition 
the persistence time, which is the time interval between two successive updates of wealth of an agent 
has been observed to have a non-trivial power law distribution. An opposite version of the model has 
also been studied where the agent with maximal wealth is selected instead of the one with minimal 
wealth, which however, exhibits similar behavior as the Minimal Wealth model. 

PACS numbers: 89.65.Gh, 87.23.Ge, 05.65.-l-b, 64.60.Ht. 



1. INTRODUCTION 

Study of the probability distribution of wealth of the 
people in a society goes back to 1897 when Pareto ob- 
served empirically that the distribution is characterized 
by a power law tail. Probability that an individual mem- 
ber has wealth more than w is given by P{w) ^ w^^" 
with a = 5/2 [l[. This type of distribution is known 
as the Pareto distribution This observation reflects 
the inherent inequality in the economic structure of the 
society. A large number of individuals are economically 
poor. In comparison the number of wealthier people is 
less and their number gradually decreases as their wealth 
increases. The cut-ofF of the distribution corresponds to 
with few very rich individuals. Consequently a sizable 
fraction of the society's net wealth is infact possessed by 
a class of top rich people consisting only a few percent 
of society's entire population. Since measuring wealth 
is difficult, in recent years distributions of income and 
tax return amounts have been studied in different coun- 
tries. For example tax return amount distribution in US 
and Japan shows a log-normal distribution in the middle 
range followed by a power law for high income people W\ , 
UK data of income shows an exponential decay which is 
followed by a power law in the high income range 0] and 
income data in Brazil for 2004 shows an almost Gaussian 
law for the low and middle income groups where as high 
salary groups are described approximately by power law 

i- 

A good amount of research effort has been devoted in 
recent times to study the wealth distribution using the 
ideas of Statistical Physics. In this description an indi- 
vidual member is called an agent. The microstate spec- 
ified by the precise description of wealth of every agent 
changes after each transaction. Their wealth change due 
to interaction among themselves. This interaction is the 



mutual bipartite trade among different pairs of agents. 
Thus the wealth distribution evolves due to such repeated 
interactions and finally assumes a time independent form. 

In a simple model the total wealth of the entire society 
has been considered as strictly conserved. In a transac- 
tion the net amount of wealth of the pair of agents is ran- 
domly re-shuffled between them and therefore the pair- 
wise interactions also maintain conservation of wealth 
0, 0] ■ Starting from an arbitrary initial distribution of 
wealth with a fixed average value the system attains a 
steady state. The steady state wealth distribution has 
been found to have an exponentially^ decaying tail. Such 
distribution was also observed in [8|. In a subsequent 
modification of the model saving propensity factor A has 
been introduced where each trader invests only 1 — A 
fraction of his wealth to the bipartite trade j9| . With a 
fixed value of A the steady state wealth distribution be- 
comes very similar to the Gamma distribution l^, 11|. 
This model was further extended by assigning quenched 
saving propensities Xi specific to each individual agent 
^12-.14] . Here one gets a power law tail with a « 2.0 
only when wealth distributions are averaged over differ- 
ent {Xi} sets 12, 1^. All these results exhibit economic 
inequality in the society. A review of all these models 
has been published in 16 1. 

Apart from these a different idea of extremal dynamics 
for the evolution of wealth distribution has been stud- 
ied by Pianegonda et. al. on an one dimensional (Id) 



lattice [17|, [18|. It is assumed that always the poorest 



agent of the society initiates a trade since he feels the 
strongest urge to raise its economic status. The trade is 
implemented by locating the poorest agent and refresh- 
ing his wealth randomly. To ensure that the model re- 
mains strictly conserved, the amount of wealth gained by 
the minimal site has been taken out equally from its two 
neighbors. Consequently this model allows an agent to 
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FIG. 1: (Color online) Time evolution of the wealth Wi at different sites of an Id lattice of size A'^ = 512. The series of snapshots 
are taken at times t = (a) 0, (b) 500, (c) 1000, (d) 500000, (e) 5000000 and (f) 15000000. The position of the poverty line at 
Wc°''-^ (512) = 0.8070 is shown by the red line. With time the wealth values gradually move up beyond the poverty line and 
stays there in the stationary state during further evolution. 



possess a negative wealth. In the stationary state of large 
systems P{w) jumps from zero to a maximal value at a 
critical value of Wc and then it decays as w increases as 
the Boltzmann-Gibbs exponential function. 

The Pianegonda model is very similar to the non- 
conservative self-organized extremal model for the eco- 
logical evolution of interacting species introduced by 
Bak and Sneppen (BS) [l^. The phenomenon of Self- 
organized Criticality (SOC) is the spontaneous emer- 
gence of fluctuations of all length and time scales in a 
slowly driven system. This concept was first introduced 
to describe the formation of a sandpile of a fixed shape 
(20| . Later the idea of SOC has been applied to a large 
number of different physical systems 2l|. In a stochas- 
tic version of the sandpile model grains are distributed 
to randomly selected neighboring sites In the SOC 
models fluctuations are described in terms of avalanches 
of activities and their size distributions assume power law 
decaying functions for large system sizes. In BS model an 
entire species is represented by a single fitness variable. 
A set of species is represented by the nodes of a graph. 
Using the spirit of Darwinian principle in each mutation 
the fitness of the node with globally minimal value is 
searched and is refreshed by a new random value. Effect 
of this propagates to few neighboring nodes which are 
also refreshed. The system eventually reaches a station- 
ary state when the fitness distribution assumes a time 
independent step like form. 

An absorbing state phase transition in presence of a 



conserved continuous local field has been studied recently 
psj . In this model a pair of sites is said to be active if 
at least one of them has energy more than a pre-assigned 
threshold. An active pair re-shuffles its net energy be- 
tween them keeping the energy strictly conserved. Be- 
yond a critical value of the threshold the number of ac- 
tive pairs fluctuates with time in the stationary state and 
the time averaged density of active pairs has been con- 
sidered as the order parameter for the problem. As the 
threshold value is tuned a continuous phase transition is 
observed from an absorbing phase to an active phase 23 1 . 
It was claimed that the critical exponents of this model 
are different from the Directed Percolation model. 

In a related model defined for the wealth distribution 
at least one of the agents in a bipartite trading is se- 
lected from a subset of agents [2J| . This subset is formed 
by agents who have wealth less than certain upper cutoff, 
the other agent being selected randomly from the neigh- 
bors of the first agent. In a transaction the net wealth of 
the pair of agents is randomly re-shuffled between them. 
The order parameter has been defined as the fraction 
of agents having wealth below a certain threshold value 
and it is claimed that the system undergoes a continu- 
ous phase transition at a critical value of the threshold 
wealth. A number of critical exponents have been mea- 
sured to characterize the transition and some of them 
are found to be close to corresponding ex pon ents in the 
Manna model of Self-organized Critically [2^ . 

In this paper we present a detailed numerical analysis 
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FIG. 2: (Color online) The system is relaxing from initial 
state to the stationary state for A'' — 2^" (black), 2^^ (red) 
and 2^^ (blue), (a) The relaxation time T has been plotted 
with the corresponding position of the poverty line wdT, N). 
(b) The same data has been plotted with deviation Wc{N) — 
Wc{T,N) on a log — log scale. The slopes are -2.30, -2.64 and 
-2.76 respectively, which on extrapolation give the value of 
the dynamical exponent z ~ 2.84. 



of the modified Pianegonda et. al. model whose trad- 
ing rule has been replaced by the stochastic re-shuffling 
of the net wealth of a pair of agents. We consider this 
model as one of the few examples of non-dissipative SOC 
systems where the entire wealth of the society is strictly 
conserved, for example the fixed energy sandpile [1^. Es- 
timation of a number of critical exponents of the modified 
Pianegonda et. al. model suggests that this model does 
not belong to the universality class of either BS model or 
Manna model. We believe that it belongs to a new uni- 
versality class perhaps because of strict conservation of 
wealth is maintained in its dynamical rules. The model 
studied by us in this paper and that in |24| are essen- 
tially same but studied from two different approaches. 
An approach which is very typical of a SOC system has 
been followed by us where the critical poverty line spon- 
taneously evolves without any fine tuning. In contrast in 
[2^ the position of the poverty line is tuned by hand to 
arrive at the critical state. 

In section 2 we describe the Minimal Wealth model 
where the trader with minimal wealth initiates the trade. 
In subsection 2.1 we discuss the relaxation dynamics of 
the system on its way to the stationary state. In sub- 
section 2.2 the correlation that evolves in the system in 
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FIG. 3: (Color online) Estimating the relaxation times, (a) 
The average wealth square {w'^{t,N)) per agent has been 
plotted with time t for N = 2* (black), 2^° (red) and 
2^^ (blue) and after a long time they approach the sta- 
tionary state value {w'^ {00, N)) . (b) The scaled deviation 
[{w^{t,N)) - {w'^{oo,N)}]N°-^^ has been plotted with the 
scaled variable tN~^'^. The relaxation exponent z — 2.70. 

the stationary state is studied. Wealth distribution in 
the stationary state is described in subsection 2.3. The 
statistics of the avalanche life-time distributions has been 
studied in section 2.4. This section ends with the study of 
persistence time distribution of individual agent's wealth 
in subsection 2.5. In section 3 we have described the 
Maximal Wealth model, which is opposite to the Minimal 
Wealth model, where the trader with maximal wealth ini- 
tiates the trade. Finally we conclude in section 4. 

2. THE MINIMAL WEALTH MODEL 

In this paper we have considered a model with a con- 
servative extremal dynamics for studying the evolution 
of wealth distribution in a society. In a bipartite trans- 
action one agent is necessarily selected as the one with 
the globally minimal value of wealth Wmin ■ The second 
agent is chosen randomly with uniform probability from 
neighbors of the first agent. This neighbor list has been 
defined in different ways for different graphs. We have 
studied this model on four different graphs, namely, (i) 
Id regular lattice with periodic boundary condition (ii) 
two dimensional {2d) square lattice with periodic bound- 
ary condition (iii) the Barabasi - Albert (BA) scale-free 
graph ^] and on an (iv) A^-clique graph. On every graph 
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FIG. 4: An exhibition of the correlation that sets in the sys- 
tem. In a single run the lattice sites imin with globally min- 
imal wealth in successive time steps are marked for a system 
of size A'^ = 512 in Id. The gradually increasing duration of 
correlation has been exhibited by time windows of increasing 
lengths: (a) 500, (b) 3000 and (c) 6000. 

the nodes represent agents and the nearest neighbors of 
each node connected by direct links constitute the neigh- 
bor hst of every agent. We report elaborately the results 
of our model on an Id lattice and mention the key results 
of the same model studied on other graphs in tables. 

The dynamics starts with N agents, each having an 
amount of wealth Wi, {i = 1, N} drawn from a uniform 
distribution with the average {w{N)) =1 irrespective of 
the system size N. The discrete time t is the number 
of bipartite transactions. At an arbitrary time t first 
the site i = imin is searched out which has the minimal 
wealth Wmin ■ The other agent j is selected randomly with 
uniform probability from the neighbors of imin- Both 
agents invest their entire amount of wealth. Therefore 
the total invested amount by both the traders is: 5ij{t) — 
Wi{t) + Wj{t). This amount is then randomly divided into 
two parts and received by them also randomly: 

w,it + 1) = e(t)%(t) Wj{t + 1) = (1) 

where e{t) is a freshly generated random fraction and 
e = 1 — e. These transactions are repeated ad infinitum. 
After some relaxation time the system reaches a station- 
ary state when the wealth distribution assumes a time 
independent form. 

In Id a linear chain of N sites with periodic boundary 
condition has been considered where the neighbor list of 
an arbitrary site i consists of the two nearest neighboring 
sites at i ± 1. Therefore the second agent j is selected 
randomly with equal probability from this list. If Wmin is 
very small then the probability that it will be replaced by 
an even smaller wealth after trade is also small. However 
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FIG. 5: (Color online) (a) The probability P{£, N) that sites 
with minimal wealth at successive time steps are separated 
by a distance £ has been plotted with £ for the system sizes 
N ^2"^ (black), 2^° (red) and 2" (blue). The slopes are -2.64, 
-2.81 and -2.89. (b) The finite-size scaling of P{£,N)N'" vs. 
IN~''^ with 77^ = 2.89 and = 1 and therefore tt = ri^/(^^ = 
2.89(5). 



this probability gradually increases as Wmin increases. As 
the sites with the minimal values of wealth are systemat- 
ically replaced, very soon all nodes with small w values 
are replaced by larger values of w resulting a vacancy in 
the small w region. This is explained pictorially in Fig. 
1. On an Id lattice with N = 512 we plot the lattice po- 
sitions i along the abscissa and the corresponding wealth 
Wi along the ordinate. As time increases a vacant region 
gradually forms for small values of w for all sites. If on 
the average the wealth of none of the agent is below a 
certain threshold value w, it is called the 'poverty line'. 
In Fig. 1 the poverty line gradually moves up with time 
and finally settles at a critical value w^"'-^ (512) = 0.8070 
at the stationary state. This behavior is the same for 
all system sizes N but with different values of Wc{N). 
Unlike the model in [iJl here the critical poverty line is 
spontaneously determined by the dynamical evolution of 
the system where no fine tuning is necessary which is the 
distinctive signature of self-organization and we will see 
in the following that the model exhibits critical behavior 
as well. 

To find the agent with minimal wealth a brute-force 
search takes cpu '-^ A^. A much faster algorithm to search 
for the globally extremal (minimal or maximal) site was 
proposed by Grassberger [27, 28] which stores the data 
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FIG. 6: (Color online) For Minimal Wealth model in Id. (a) 
The wealth distribution P{w, N) in the stationary state for 
iV = 2'^ (black), 2" (red) and 2^^ (blue), (b) A finite-size 
scaling using P{w, N)N~°-°'^ vs. [w-Wc{N)]N°-'^^ for system 
sizes iV = 2" (black), 2^^ (red) and 2" (blue). 

in a Hierarchical structure. This takes cpu ^ In TV. We 
have used this method for Id, 2d and for A^-clique graphs. 
For BA graphs we used the brute-force algorithm. 

2.1 Relaxation to Stationary State 

We first estimate the relaxation time required for the 
system to reach the stationary state. During this relax- 
ation period the wealth distribution gradually changes 
starting from its initial uniform distribution to its time 
independent form in the stationary state. The relaxation 
time has been estimated as a function of deviation of the 
poverty line from its critical value in two ways. 

For a given system size N we have simulated our model 
up to 10^ time steps and calculated the wealth distribu- 
tion P{T,w, N) at 100 time instants T at the interval of 
lO"' steps. These distributions are calculated for a sample 
size of ~ 10^ independent runs. The value of the poverty 
line Wc (T, N) at time T is determined by the value of 
w where P{T,w, N) is the maximum. This estimation 
is done for all values of T. In Fig. 2(a) we plot on a 
lin - lin scale T vs. WciT,N) for N = 2^°, 2" and 
2^^. We see that in each case the relaxation time T di- 
verges as WciTjN) approaches its stationary state value 
Wc{N). These data have been replotted in Fig. 2(b) us- 
ing a log — log scale as T vs. Wc{N) — Wc{T, N) where we 




FIG. 7: (Color online) Gaussian fit of the wealth distribution 
P{w, N) in the stationary state for seven different system sizes 
iV = 2^^ to 2" in Id. The sequence starts with the distribution 
for A'' = 2^ at the extreme left and is shifted to the right by 
0.2 when system size is multiplied by a factor of 2. 

have used Wc{N) = 0.8242, 0.8375 and 0.8383 to obtain 
the best straight line plots. The slopes of these straight 
lines are 2.30, 2.64 and 2.76 respectively which are then 
extrapolated with ]\[-^-^03 obtain the exponent z as 
T{N) cx [wc{N) - Wc{T, N)]~'- with z « 2.84. 

In a second method we calculated {w'^{t, N)) with time 
t starting from its initial value when the distribution is 
uniform. After a long time this quantity saturates to its 
stationary value {w'^{oo, N)). In Fig. 3(a) we show the 
plots of {'w'^{t, N)) vs. i on a log — log scale. In Fig. 3(b) 
[{w^{t, N)) - (w2(oo, N))]N°-^^ has been plotted with the 
scaled value of time tN~'^-^ using (w^(oo,iV)) — 0.057, 
0.0607 and 0.061 for N = 2^ 2^° and 2^2 respectively. A 
nice data collapse has been obtained with the following 
scaling form 

[{w^{t,N)) - (ii;2(oo,iV))]7V0-55 ^ J-,(tiV-2-5) (2) 

where the scaling function J^z{x) varies as x~^^^ for small 
X. A direct measurement of the slope of the scaled plot 
gives an estimate of 1/z = 0.37 which corresponds to 
z = 2.70. We conclude an average value of z — 2.77(7). 

2.2 Correlation in the Stationary State 

Starting from an uncorrelated wealth distribution the 
system becomes more and more correlated as time passes. 
This is reflected in the fact that the probability of occur- 
rence of the minimal sites close to each other in successive 
time steps gradually increases. At the early uncorrelated 
stage the position of the minimal site at the next time 
step is likely to be anywhere in the lattice. However as 
time increases, the poverty line moves up, consequently 
Wmin increases and the probability that the minimal site 
at the next time step occurring at the same site or at the 
neighboring updated site also increases. This is shown 
in Fig. 4 using a Id lattice of iV = 512. For a single 
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FIG. 8: (Color online) Estimation of the critical poverty 
line U'c(oo). The w^"'''^ {N) are the upper plot (in red) and 
w^'°^''{N) are the lower plot (in blue). When they are extrap- 
olated with N'" with K.{hal}) = 0.649 and K,{slope) = 0.657, 
it)c(oo) is obtained as 0.8174 and 0.8176 respectively. We con- 
clude Wc(oo) = 0.8175(2). 



run it is observed that the locations of Wmin are quite 
random (Fig. 4(a)). However as time evolves these posi- 
tions gradually become more correlated (Figs. 4(b) and 
4(c)). In general one can consider the successive jumps 
of imin positions constituting a Levy flight random walk 
[jgj . We see below that indeed their flight lengths follow 
a power law distribution. 

The correlation in the stationary state is quantita- 
tively measured by the probability distribution P{1) of 
the distance of separation ^ between successive minimal 
sites using periodic boundary condition. This distribu- 
tion measured in the stationary state has been plotted 
in Fig. 5(a) for different system sizes A'' = 2^,2^" and 
2^^. The value of P{t} for £ = and 1 are approximately 
0.4575(1) and 0.4820(1) and then it decreases as a power 
law P{t) ^ with increasing £. A direct measurement 
of slope gives tt « 2.89. Fig. 5(b) exhibits the finite-size 
scaling of the same data when the £ and P{£, N) axes are 
scaled as: 



(3) 



where J- nix) is a universal scaling function with the scal- 
ing exponents rj^r — 2.89 and Ctt = 1- From this scaling 
analysis we get n = t^^/Ctt = 2.89(5). 

There exists a spatial correlation too in this model. 
A two point correlation function has been measured in 
the stationary state. The average correlation between 
two sites situated at a distance of separation x has been 
defined as: 



C{x) = {w{0)w{x)) ~ (w) 



(4) 



where (w) is always set equal to unity. We assume a 
power law decay for the correlation, i.e., C{x) ~ for 
X — >■ oo. For Id a plot of (not shown) C{x) vs. a; on a 
log — log scale indicates a power law for large x values. 
However considerable variation of slopes exists for system 
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FIG. 9: A portion of the time series of minimal wealth values 
in successive time steps is shown for a Id system with TV — 256 
and Wc(256) = «)c'''^(256) = 0.8010 in the stationary state. 
It shows that depending on the value of Wo an avalanche can 
be broken into a hierarchy of avalanches. For this run the 
system has been rela^ced for the initial tskip — 5 x 10* time 
steps. 



sizes TV = 2*,2i" and 2^2. The slopes are: -1.17, -1.29. 
-1.34 respectively which extrapolates to a value of x = 
1.5(1) in the limit of — > oo. Our estimate for x in 2(i 
is 2.2(2). 



2.3 Wealth Distribution in the Stationary State 

Next we estimated the probability density distribution 
P{w, N) of wealth in the stationary state of the system of 
size N . The distribution grows very rapidly near wdN) 
for all values of N and then decays very fast. In Fig. 6(a) 
we show the plot of P{w, N) vs. w tor N — 2^, 2^° and 
2^^ . All of them have similar variations but with increas- 
ing system size the curves gradually become sharper. 
Therefore we tried a finite-size scaling analysis in Fig. 
6(b) for the growing region and for N = 2^^, 2^^ and 2^^. 
A nice data collapse is observed when axes are scaled and 
P{w, iV)iV-0-0'* has been plotted with [w - WciN)]N°-''^ . 

The functional form of the decay of the probability dis- 
tribution P(w, N) has been studied right after the max- 
imal jump at WciN). This part fits very well with the 
Gaussian form: 



P{w,N) 



A{N) 



2T:a{N) 



exp[ 



2a^{N) 



(5) 



In Fig. 7 we showed P{w,N) vs. w on a lin — lin scale 
for seven different system sizes: TV = 2^, 2®, 2^'^. In 
each case the fitting curve is indistinguishable from the 
data. We observe a systematic variation of A{N), ii{N) 
and (t(N) with system size N. For example A{N) w 
511.5 - 1828iV-o-294^ ^(7v) « -1.10 -t- 62.5iV-i-2i5 ^^^^ 
a{N) w 0.64 + 0.4297V-0-436. 

The precise value of Wc{oo) is calculated by extrapo- 
lating Wc{N) values which are calculated using the fol- 
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FIG. 10: (Color online) For Id Minimal Wealth model, (a) 
The avalanche size distribution for A*' = 2^^ (black), 2^^ (red) 
and 2" (blue), (b) A finite -size scahng of the data in (a) 
with scaling exponents r^r = 2.44 and (^t — 2.18 giving the 
avalanche size exponent r — rj-r/^r ~ 1.12(1). 

lowing two methods. We have seen in Fig. 6(a) that 
the probability distribution of P{w,N) becomes increas- 
ingly steeper with increasing N. For a certain size N we 
defined w'^"'^^ {N) as the value of w for which P{w,N) 
is half of its maximum value. In a second method the 
Wc{N) has been calculated in the following way. A pair 
of successive points on the P{w,N) vs. w curve which 
has the largest slope is found out. A straight line join- 
ing these two points is then extrapolated to meet the w 
axis at wf°P''{N). The pair of values of w'^^^f [N] and 
wf°P''{N) for eight N values 2^ to 2^^ arc then extrap- 
olated with N''^. A least square fit of straight line has 
been done for trial values of k starting from 0.20 to 1.20 
at an interval of 0.001 and the errors have been calcu- 
lated. The errors are minimal for K{half) = 0.649 and 
K{slope) = 0.657. Using these two values of the expo- 
nent n we extrapolate Wc{N) values with iV"" in Fig. 8 
to meet the Wc{N) axis at 0.8174 and 0.8176 respectively. 
We conclude a value for Wc{oo) for Id model as 0.8175(2). 

2.4 Avalanche Size Distribution 

In the stationary state successive values of minimal 
wealth Wmin fluctuates with time. If a certain reference 
wealth is fixed by hand at w = Wo then the successive 
Wmin appear below and above Wo hne. One defines a Wq- 
avalanche as the sequence of successive bipartite trades 
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FIG. 11: The average size of the avalanches {s{wc{N), N)) 
right at the poverty hne Wc{N) has been plotted with system 
size N on a. log — log scale. The slope of this curve gives the 
exponent /3 = 1.92(2). 

whose Wjnin values are smaller than Wq. The size s of the 
avalanche measures the duration of the avalanche i.e., at 
times t and t + s + 1 the Wmin > Wo, whereas at every 
time step from t + \ to t + s the Wmm < Wo- When 
Wo is set equal to Wc{N) it is called a critical avalanche. 
This is explained in Fig. 9 where part of the Wmin time 
series for N = 256 and Wc(256) = u;J°'/(256) = 0.8010 
is displayed discarding the initial tskip — 5 x 10^ time 
steps. For Wo — 0.71 an avalanche of size 53 breaks into 
two avalanches of sizes 25 and and 27 when Wo is reduced 
to 0.66. On further reduction of Wo to 0.56 these two 
avalanche break into even smaller avalanches of sizes 18, 
6 and 10, 16 respectively. Thus any avalanche can be 
splitted into a hierarchy of smaller avalanches if Wo value 
is lowered [s^. On the other hand if Wo is raised the 
average avalanche size increases and becomes infinite at 
certain value of Wo- 

At the critical point the distribution of the avalanche 
life-times has a power law tail in the limit of A'^ — ?> oo: 
P(s, oo) ^ s^'^. In the stationary state we used Wo — 
^ha;/j-jY) measured life-times of a large number of 
avalanches for different system sizes and plot the proba- 
bility distributions P{s, N) vs. s using a log — log scale 
in Fig. 10(a). Each curve has a straight portion in the in- 
termediate regime of the avalanche sizes and this regime 
becomes gradually larger on increasing A^. The direct 
measurement of slopes in the scaling region gives t{N) 
= 1.086, 1.091 and 1.096 for A^ = 2^^ 2'^^ and 2^4 respec- 
tively. A finite-size scaling is very much suitable with the 
following scaling form: 

P(s, A)A^''- cx Tr{sN-^-) (6) 

where the scaling function J-t{x) ^ x^'^ in the limit of 
X — and J>(a;) approaches zero very fast for a; >> 1. 
The exponents r]r and fully characterize the scaling 
of P{s, N) in this case. An immediate way to check the 
validity of this equation is to attempt a data collapse by 
plotting P{s,N)N'^^ vs. s/N^^ with trial values of the 
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lo' 10^ lo' lo" lo' 10^ lo' 




FIG. 12: (Color online) (a) The average value of the avalanche 
size {s{wo,N)) has been plotted with deviation Wc{N) — Wo 
from the poverty line Wc{N). The system sizes used are 
iV = 2* (black), 2" (blue), 2^2 (green) and 2" (red). The 
value of 7 obtained by extrapolation of slopes is 2.67. (b) Fi- 
nite size scaling analysis of the data in (a) is shown. The 
scaling exponents rj-y — j3 = 1.92 and (.y = 0.73 gives 
7 ~ ViKi ~ 2.63. (c) Plot of the data in (a) but with 
Wc{N) — w^"^ {N). The slopes for four different system 
sizes in (a) on extrapolation gives 7 = 2.66. We conclude 
7 = 2.65(5). 



scaling exponents. For Id the values for obtaining the 
best data collapse are found to be rjr = 2.44 and Ct — 
2.18 (Fig. 10(b)). The life-time distribution exponent 
for Id is therefore r — titIC^-t ~ 1.12(1). 

Next we calculated the average value of avalanche life- 
times {s{wc,N)) right at the critical poverty line. In 
Fig. 11 we plot this quantity with system size iV on a 
log — log scale. The plot fits excellent to a straight line 
and its slope gives the value of the exponent /3 w 1.92(2) 
m: {s{wc,N)) ~ 7V^. Ass uming the distribution P(s, N) 
of avalanche sizes holds good up to a cut-off Smax ~ N'^^ 
one gets a scahng relation /? ~ Cr (2— r) and our estimates 
of /3 = 1.92, = 2.18 and r = 1.12 satisfy this relation 
very closely. 




FIG. 13; (Color online) The probability distribution P{tp, N) 
of persistence times tp at the stationary state, (a) The finite- 
size scaling of the distribution in Id for A'^ = 2^" (black), 2^^ 
(red) and 2^* (blue). From the scaling exponents rje = 3.41 
and (e = 2.21 the persistence exponent 6 = rje/C.e ~ 1.543 is 
obtained, (b) The finite-size scaling of the distribution in 2d 
for L = 2* (black), 2® (red) and 2^° (blue). From the scaling 
exponents rje = 3.49 and C,e = 2.74 the persistence exponent 
9 = rje/Ce ~ 1.274 is obtained. 

The size of the Wo-avalanches are smaller when Wo < 
WciN) and we have studied how the average avalanche 
size grows as the deviation {wc{N) ~ Wo) decreases. Sim- 
ilar to the BS model we assume {s(wo)) ~ [wc — Wo\~'^ 
for iV — ^ 00. We measured the average size {s{wo, N)) of 
the Wo-avalanches for different system sizes N and plot- 
ted them with Wc{N) — Wo in Fig. 12(a) with wdN) = 
^haif^j\j-y all plots on log — log scale the curves are 
horizontal as deviation Wc{N) — w is very small. How- 
ever as the deviation increases the curves become straight 
with negative slopes -1.98, -2.15, -2.28 and -2.38 for for 
N — 2^, 2^°, 2^^ and 2^"* respectively. These values when 
extrapolated with 7V~°-^°® give 7 — 2.67 for iV — >• 00. 
Again a finite-size scaling has been possible as shown in 
Fig. 12(b): 

{siwo,N))N-'^-' cxT^i[wciN)~Wo]N^-) (7) 

where J^y(a;) is another scaling function. From this data 
collapse the scaling exponents t]^ = j3 ^ 1.92 and — 
0.73 with 7 — i]-ylQ'y ~ 2.63 is obtained. 

For every system size N there is a value of Wo = 
w™^ {N) so that when Wo is raised to this value the 
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FIG. 14: (Color online) Values of wealth Wi of different agents 
in the stationary state are plotted with their positions i along 
a Id lattice for the (a) Maximal Wealth model, the red line is 
at w^"''^ (512) = 1.3924. (b) mixture of the Minimal Wealth 
and Maximal Wealth models with probability p — 1/2, the 
red lines correspond to peaks at ?j;c(512) = 0.7559 and at 
1.3574. 



avalanche size becomes infinite. This implies that if we 
plot the data in Fig. 12(a) with respect to w"''^ [N) — Wo 
then we should be able to see the divergence of aver- 
age avalanche size instead of saturation of the avalanche 
sizes. We plot this in Fig. 12(c) using wl"f (N) = 0.8167, 
0.8169, 0.8170, 0.8172 for N = 2^, 2^", 2^2 and 2^'^ re- 
spectively again on a log — log scale. Each curve is a 
straight line but with different slopes: -2.31, -2.43, -2.51 
and -2.56 respectively. When these slopes are extrapo- 
lated with iV~°-^^ the extrapolated value for — >■ oo 
is -2.66. Our conclusion for the value of the exponent 
7 = 2.65(5). 



2.5 Persistence of Wealth in the Stationary State 

The time interval between two successive updates of 
wealth of an agent is known as the persistence time tp. 
Different agents have to wait for different amounts of 
times in general. More specifically agents having small 
amount of wealth have to wait for very little times. On 
the other hand potentially rich agents have to wait long 
enough times. In the stationary state we measure the 
persistence times for all sites of the lattice and use this 
data to plot their probability distribution. More precisely 
we set a clock to each site. Whenever there is a change 



of wealth at this site the time is noted and the clock time 
is reset to zero. At the stationary state we collect a large 
number of persistence time data and use these data to 
calculate the persistence time distribution. 

We assume a power law variation of the persistence 
time distribution as P{tp) ^ V 



For finite size systems the distributions P{tp,N) vs. i,p 



in the limit of — > cxi. 

are plotted on a log — log scale (not shown) and the direct 
measurement of slopes give the 0{N) values for finite size 
systems. These values are extrapolated as N'-^-'iQi 
obtain 9 ~ 1.539 for — > oo in Id. In a similar analysis 
for a 2d square lattice of size L using an extrapolation 



with respect to L 



-0.565 



we get 



1.25 for L — >• oo. 



Persistence exponents are also obtained by the finite- 
size scaling analysis. In Fig. 13(a) we show the scal- 
ing plot of P(tp,N)N'^o vs. tpN-"^" with rie = 3.41 and 
Ce = 2.21. This gives 9 = T]e/Ce = 1-543 in Id. Similar 
scaling analysis in terms of the system size L in 2d square 
lattice has been performed with ije = 3.49 and — 2.74 
which gives 9 = 1.274 (Fig. 13(b)). Averaging 9 values 
obtained from direct measurement and scaling analysis 
we conclude 9 = 1.541(10) in Id and 9 ^ 1.262(10) in 2d. 



3. THE MAXIMAL WEALTH MODEL 

Next we studied the Maximal Wealth model where one 
agent is necessarily the agent with maximal wealth. The 
other agent being selected randomly with uniform prob- 
ability from the neighbors of the first agent. Random 
re-shufSing of wealth takes place in the same way as in 
the Minimal Wealth model. 

In Fig. 14(a) we plot again for the Maximal Wealth 
model the values of wealth Wi of different agents at a 
certain instant of time in the stationary state with their 
positions i along an Id lattice of size N = 512. In contrast 
to the similar plot of the Minimal Wealth model in Fig. 
1 here an upper cut-off for the wealth has been visible at 



,,half 



(512) = 1.3924. 



In this case the stationary state wealth distribution 
P{w,N) takes an opposite shape (Fig. 15(a)). A crit- 
ical wealth Wc{N) exists here as well. P{w,N) takes 
a Gaussian form elevated by a constant term c{N) for 
w < Wc{N), whereas for w > Wc{N) it sharply de- 
creases to zero. The parameters of the Gaussian func- 
tion (Eqn. (4)) are different for different N and they 
vary very systematically with N as: A{N) w 30.48 — 
133.357V-0-36i^ ^(7V) ~ 3.74 _ 6MN-°-^^ and a{N) w 
1.047 + 0.5727V-''-32i g^^^ ^Yie constant c{N) w 0.031 + 
4.265A^~"-^3_ 

The critical value of wealth Wc{oo) in the asymptotic 
limit of system sizes has been estimated by the same 
method as used for the Minimal Wealth model. The 
w^°'^f {N) and wf°P^{N) values have been calculated for 
N = 2^, 2^°, 2^2 and 2", extrapolated with N^o.bse ^^^^ 
^-0.620 g^^^ i-j^g asymptotic values are 1.3610 and 1.3608 
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Minimal Wealth model 


BS model 


Manna model 


Id 


2d 


BA graph 


TV-clique 


Id 


2d 


Id 


2d 




0.8175(2) 


0.6887(2) 


0.6444(2) 


0.6076(2) 


0.66702(8) [271 


0.328855(4) [30] 


0.89199(5) [31] 


0.68333(3) [31] 


T 


1.12(1) 


1.29(1) 


1.50(1) 


1.50(1) 


1.073(3) [22 


1.245(10) \30i 


1.112(6) [32] 


1.273(2) [32| 


7 


2.65(5) 


1.58(5) 


1.02(5) 


1.00(5) 


2.70(1) [30] 


1.70(1) [30| 








2.89(5) 


3.94(5) 






3.23(2) [30] 








Z 


2.77(7) 
















P 


1.92(2) 


0.95(2) 


0.52(2) 


0.50(2) 






2 


2 


e 


1.541(10) 


1.262(10) 


1.00(1) 


1.00(1) 











TABLE I: Values of different exponents of Minimal Wealth model are compared with those of existing models in the literature. 
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FIG. 15: (Color online) (a) The wealth distribution P{w, N) 
vs. w for the Maximal Weahh model for N = 2^ (black), 2^° 
(red), 2^^ (green) and for 2^'* (blue), (b) Finite-size scaling of 
the avalanche size distribution P{s, N) at the critical thresh- 
old Wc{N) for N ^2^ (black), 2^° (red) and 2^^ (blue). The 
scaling exponents are rjr = 2.39 and = 2.10 which gives 
the exponent r = rir/C,T ~ 1.14(1). 



FIG. 16: (Color online) (a) Mixture of Minimal Weahh and 
Maximal Wealth models with probabilities p and 1 — p re- 
spectively for iV = 2^ and for p = 0.02 (black), 0.08 (red), 
0.2 (green), 0.4 (blue), 0.6 (magenta) and 0.9 (maroon), (b) 
Variation of the ratio r(p) of heights of the right peak and the 
left peak with probability p. 



respectively. We conclude Wc(oo) = 1.3609(2) for Id. A 
snnilar analysis gives Wcioo) = 1.7076(2) for 2d square 
lattice, 1.8895(2) for the BA graph and 1.9998(2) for the 
TV-clique. 

It may appear that the Minimal Wealth and Maximal 
Wealth models should be symmetric about the average 
wealth per agent which we have set at {w) = 1. We have 
seen above that this indeed not the case since Wc values 
are 0.8175 and 1.3609 for the Minimal Wealth and Maxi- 
mal Wealth models respectively. The symmetry between 
these two models are broken by the presence of a rigid 
wall &% w — Q which means that negative value of wealth 
of an agent is not allowed. 

The avalanche size distributions have been studied as 



well. A finite-size scaling of these distributions has been 
done and are plotted in Fig. 15(b) using log — log scale 
as before for N = 2^, 2^° and 2^^. The scaling exponents 
are rjr = 2.39 and Ct = 2.10 respectively giving the value 
of the avalanche size exponent r = ryr/Cr ~ 1.14(1). 

Finally we have studied a mixture of the Minimal 
Wealth and Maximal Wealth models. At every bipar- 
tite trade the first agent with minimal wealth is selected 
with probability p or with maximal wealth with prob- 
ability 1 — p. The second agent is selected with uni- 
form probability from the neighbors of the first agent. 
A snapshot of the individual wealth for p = 1/2 at the 
stationary state for different agents has been shown in 
Fig. 14(b). Here the wealth values are restriced within 
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a 'wealth-band' with sharp cut-offs at a high and a low 
end. Consequently the shape of the wealth distribution 
P{w,p) has double peaks for all N and we plot the distri- 
bution in Fig. 16(a) for N = 2^. Portion of the distribu- 
tion between the peaks fits excellent to elevated Gaussian 
distributions with different parameter values for different 
values of p in the range between 1/2 and 1. On the two 
sides of this region the distribution decays to zero very 
fast. The figure shows the plot of P{w,p) vs. w for p = 
0.02, 0.08, 0.2, 0.4, 0.6 and 0.9. In Fig. 16(b) we plot the 
ratio r{p) of the heights of right peak and the left peak 
with the probability p. 

4. CONCLUSION 

Social inequality in terms of economic strengths is 
ubiquitous for the people of all countries. Perhaps this 
inequality acts as the major driving force behind the ad- 
vancement of society. Consequently the mechanism that 
establishes this inequality in a society is an important 
issue and attracted the attention of researchers over the 
last century. Here we have studied a modified version 
of the conservative self-organized extremal model intro- 
duced by Pianegonda et. al. which is motivated by the 
wealth distribution in a society. In this model the entire 
wealth of the society is strictly conserved. It evolves by 
a trade dynamics that takes the society from equality (or 
any other initial wealth distribution) to a stationary state 
where strong social inequality is present. The dynam- 
ics is an infinite sequence of stochastic bipartite tradings 
where one of the agents has the globally minimal value of 
wealth, the other one being selected randomly from the 
neighbors of the first agent. Our numerical study reveals 
that this model is one of the simplest models of Self- 
organized Criticality where the stationary state is non- 
ergodic. This model is very similar to the self-organized 
critical Bak-Sneppen model for the ecological evolution of 
interacting species. Using numerical simulation we have 
estimated a number of critical exponents for this model 
on an Id regular lattice, 2d square lattice, the Barabasi 
- Albert scale-free graph and on the A^-clique graph. We 
present evidences which suggest that this model does not 
belong to the universality class of either the Bak-Sneppen 
model or the Manna model of Self-organized Criticality. 
This model belongs to a new universality class perhaps 
because of strict conservation of wealth is maintained in 
its dynamical rules. 
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